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ABSTRACT 

Inspired by the f(R) non-linear massive gravity, we propose a new kind of modified gravity 
model, namely f(T) non-linear massive gravity, by adding the dRGT mass term reformulated in 
the vierbein formalism, to the f(T) theory. We then investigate the cosmological evolution of f(T) 
massive gravity, and constrain it by using the latest observational data. We find that it slightly 
favors a crossing of the phantom divide line from the quintessence-like phase (wde > —1) to the 
phantom-like one (Wde < — 1) as redshift decreases. 
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I. INTRODUCTION 

Einstein’s theory of general relativity (GR) has achieved great success in explaining various gravi¬ 
tational phenomena since its birth, and till today it can pass the Solar System’s tests to a very high 
precision. However, on scales larger than the Solar System, there are some astrophysical observations 
having not yet been fully understood or explained in the context of GR. Maybe the most urgent among 
them are the presence of an invisible (dark) matter component in our universe, first discovered by Zwicky 
in the 1930s [1, 2], as well as the late-time acceleration of our universe confirmed by various cosmological 
observations [3-11]. Although by extending the standard model of particle physics, there exist various 
promising candidates for dark matter (see e.g. reviews [12-14]), so far there has been no natural and 
satisfactory explanation for the cosmic acceleration by introducing some kinds of mysterious cosmologi¬ 
cal fluids with negative pressure called dark energy (DE), since they all suffer from the fine-tuning [15] 
and/or the cosmic coincidence problems [16-18]. Unless all of these problems are eventually solved on the 
particle physics side, the severity of the phenomenological problems challenges GR as the ultimate and 
complete theory of gravity. Furthermore, difficulties arise when trying to search for a quantum theory of 
gravity, since GR is not renormalizable [19, 20] and thus cannot make meaningful physical predictions. 
It is, therefore, reasonable to investigate the possibility of modifying Einstein’s gravity theory. 

Among the modified gravity theories, f(R) theory (see e.g. [21, 22, 91, 92] for recent reviews) serves 
as one of the most straightforward and popular generalizations of GR by extending the Ricci scalar R in 
the Einstein-Hilbert action to a general function f(R). It has been shown that f(R) gravity can explain 
the present cosmic acceleration without the need of DE [23-32, 93, 94]. But in general, the resulting field 
equations of motion are 4th order in f(R) gravity. This feature makes f(R) theory quite hard to analyse. 
In comparison with f{R) gravity, there is another kind of modified gravity, namely f(T) theory, which 
has been extensively explored in the literature (see e.g. [33-44]). Here T is the torsion scalar constructed 
from the curvatureless Weitzenbock connection [45]. It is well known that f(T) theory is based on the 
teleparallel space-time, in which vierbein field is used as the basic dynamical variable instead of metric 
in GR. One advantage of f(T) theory in contrast to f(R) theory is that, the field equations arising from 
f(T) theory are 2nd order, and thus simpler than f(R) case. 

On the other hand, as a classical field theory, GR propagates a non-linear self-interacting massless 
spin-2 field. Therefore, another way to modify GR is to give a tiny mass to the graviton. This class of 
theories is known as massive gravity in the literature. Actually, massive gravity has a long and elusive 
history. The first attempt to give a mass to graviton dates back to the year of 1939 when Fierz and Pauli 
proposed a unique ghost-free mass term to gravity in the linear perturbation level [46]. The linearised 
massive gravity propagates 5 degrees of freedom (DoF). And then, in 1970s, van Dam-Veltman-Zakharov 
(vDVZ) discontinuity was found. It was shown that, when the mass of graviton m —> 0, Fierz-Pauli theory 
would not converge to massless case [47, 48], thus giving different predictions than those GR. Almost at 
the same time, Vainshtein [49] showed that the linear approximation loses its validity around massive 
sources below the Vainshtein radius and non-linear effects should be considered. The vDVZ discontinuity 
can then be avoided by this well-known Vainshtein screening mechanism. However, it was claimed that 
any non-linear extensions of Fierz-Pauli theory will generally have 6 DoF, with the extra DoF being a 
Boulware-Deser (BD) ghost [50] . The construction of a well-defined non-linear theory of massive gravity 
without the BD ghost has been a challenging problem over decades. Recently, it was shown that BD 
ghost may be removed order-by-order in a careful construction of mass term up to a certain decoupling 
limit [51]. The full non-linear massive gravity was resummed by de Rham, Gabadadze and Tolley (dRGT) 
in [52] , and the absence of BD ghost was further confirmed by using Hamiltonian constraint analysis in 
[53, 54]. Although dRGT theory exists accelerating solutions [55-58], it suffers from some instabilities 
at the cosmological perturbation level [59-61]. Since then dRGT theory has been extended in various 
approaches, e.g. the mass varying massive gravity [62, 63] and the quasi-dilaton massive gravity [64-66]. 

Recently, a new extension of massive gravity, namely f(R) non-linear massive gravity, has been pro¬ 
posed [67, 68]. It is shown that this model not only shares no linear instabilities around a cosmological 
background, but can also describe the inflation and the recently cosmic acceleration in a unified frame¬ 
work [67, 68, 95]. Inspired by the f(R) non-linear massive gravity, we propose a new kind of modified 
gravity theory by adding the non-linear dRGT mass term formulated in the vierbein formalism, to the 
f(T) theory. We might call this theory f(T) non-linear massive gravity, or simply f(T) massive gravity. 
One merit of our model compared with the f(R) massive gravity is that, the field equations arising from 
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f(T) massive gravity are 2nd order, and thus much easier to analyse. Moreover, f(R) massive gravity, 
as well as dRGT theory, is generally constructed in the metric formulation. However, in the metric 
formalism, the mass term is difficult to analyse, since it contains the square root of the metric. Another 
merit of our theory is that, in teleparallel space-time, we use vierbein instead of metric as the basic 
dynamical variable. Since vierbein can be viewed as the square root of the metric in some sense, the 
vierbein formalism in our theory can greatly simplify the calculation. 

After presenting the action of our model, we then investigate the cosmological evolution of f(T) massive 
gravity. The rest of this paper is organized as follows. In Sec. II, we briefly review the f(T) gravity. 
In Sec. Ill, we add the ghost-free dRGT mass term reformulated in the vierbein formalism to the f(T) 
sector, thus giving the action of f(T) massive gravity. The evolutionary equations for the flat Friedmann- 
Robertson-Walker (FRW) cosmology are given as well. In Sec. IV, we focus on two concrete examples, 
namely the power-law and exponential types of f(T) massive gravity and investigate their cosmological 
evolution. We constrain the model parameters by the recent cosmological data and study the phantom 
crossing behavior. In Sec. V, we give some concluding remarks. 


II. f(T) GRAVITY 


/(T) gravity, which is a generalization of the teleparallel gravity originally proposed by Einstein [69-73], 
has received great interest in the literature recently. Following [33, 74, 75], we now briefly review the f(T) 
theory. f(T) gravity is built upon teleparallel space-time. In teleparallel space-time, the basic dynamical 
quantity is a vierbein field e a = e a p dp : with Latin indices a, b, ■ ■ ■ = 0,--- ,3, and i,j,-■■ = l,-- - ,3, 
Greek indices /r, zz, • • • = 0, • • • ,3, and d M coordinate bases. We also note that the Einstein summation 
notation for the indices is used throughout this paper. The vierbein is an orthonormal basis for the 
tangent space at each point x M of the manifold, namely e 0 • e;, = rj a b, with rjab = diag (—1, 1, 1, 1). The 
metric tensor can then be expressed in the dual vierbein e a M as 

SV( x ) = rj ab {x)e h v (x). (1) 

Rather than using the torsionless Levi-Civita connection in general relativity (GR), teleparallel space-time 
uses the curvatureless Weitzenbock connection v , which is defined by 


r A = P x r) P a 
L pv — c a V 


( 2 ) 


Note that the lower indices /r and v are not symmetric in general, thus the torsion tensor is non-zero in 
the teleparallel space-time. The Weitzenbock torsion tensor is defined by 


T- 


fill 


= r A - r 

- 11 fl 


fill 


= e 0 A ( d u e a ^ - dpe a u ) . 


( 3 ) 


In teleparallel gravity, the gravitational action is given by the torsion scalar instead of the the Ricci scalar 
in GR. The torsion scalar is basically the square of the Weitzenbock torsion tensor, and is defined by 


rp _ np rp p,l/ —/pp rp pi/ , ^_*pp rpl/p rpp rpl/p 

^ ^ Ills' p ^ 1 fill p ' 2 1 P HP 1 V 


with the tensor S p ^ v given by 


S p p v = \ (n* - V + r„/) + \^T° va - l -KT 


per ’ 


( 4 ) 

( 5 ) 


In f(T) gravity, the gravitational field is driven by a Lagrangian density T + /(T), with f(T) a function 
of T, and the action reads 


5 = -Ml Je[T + f(T)} d 4 x + S m [e/, Xm ], (6) 

where e = det (e a M ) = y/^g, Ml = (167rG) _1 , with g the determinant of the metric gp U and G the 
Newtonian constant. Note that we have used the units in which the speed of light c = 1, and the reduced 
Planck constant h = 1. Here, S m [e a p - Xm] is the matter part of the action, and denotes all matter 
fields collectively. If we set f(T) = 0, then the action gives an equivalent description of the space-time 
as GR. 
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III. f(T) NON-LINEAR MASSIVE GRAVITY 

In this section, we will add the dRGT mass term reformulated in the vierbein form to the f(T) sector, 
thus giving the action of /(T) massive gravity. We then derive the evolutionary equations for the flat 
FRW cosmology. 


A. The action 

Since teleparallel gravity equals to GR at the field equations level, here we give an equivalent description 
of dRGT massive theory by adding the mass term to teleparallel gravity. The mass term firstly proposed 
in dRGT theory, can be greatly simplified when transforming to the vierbein formalism. We refer to [76] 
for a detailed derivation. By adding the mass term to teleparallel gravity, one will get the teleparallel 
version of ghost-free dRGT massive gravity as 

S = —Mp f eTd^x — Mpin 2 f ^ e/3 n 5 n (E) [e/, Xm], (7) 

^ ** n—0 

where f3 n are free parameters, E denotes the matrix of vierbein ej 1 and S'n(M) are the so-called elementary 
symmetric polynomials [76]. For an arbitrary 4x4 matrix M, the first few of S n (M) can be written in 
the form [76] 


S i(M) = 

[M], 

(8) 

S 2 (M) = 

i ([M] 2 - [M 2 ]) , 

(9) 

S 3 (M) = 

^ ([M] 3 - 3[M][M 2 ] + 2[M 3 ]) , 

(10) 


with [M] the trace of the matrix M. By generalizing the torsion scalar T to a function T + f{T) in the 
action (7), we have the /(T) non-linear massive gravity as 

5 = -M 2 fe[T+ f(T)} d 4 x - Mlm 2 / ^ ep n S n (E) d 4 x + S m [e/, Xm ], (H) 

n= 1 

Note that we have absorbed the cosmological constant arising from the mass term to the f(T) sector. 
So, there are only three free parameters /3i, /3 2 , and /?3 in the mass term. We also note that when 
f(T) = 0, the action (11) is just the teleparallel gravity plus a dRGT mass term reformulated in the 
vierbein formalism, thus being equivalent to the original dRGT theory at the field equations level. So, 
if f(T) = 0, then the BD ghost will not show up in the action (11). However, if ,/(T) is a non-trivial 
function of T, then additional DoF will be introduced. Therefore, we should perform Hamiltonian analysis 
to confirm that the action (11) is ghost-free. In this paper we focus on the cosmological behavior of f(T) 
non-linear massive gravity, and the ghost problem or other instabilities that our theory may suffer from, 
will be considered in our future work. 


B. Equation of motion 

We now consider a spatially flat FRW space-time 

ds 2 = —N 2 (t) d 2 r + a 2 (r) 8ijdx l dx\ (12) 

where a(r) is the scale factor and N(t) is the lapse function. This metric arises from the diagonal dual 
vierbein 

e% = diag (N(t), a(r), a(r), a(r)) (13) 
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through Eq. (1). Then we can easily get ej 1 = diag (iV 1 , a 1 , a 1 , a x ), e = Na 3 , and hence 

s,m = 4+w. 

a z Na 


s 3 { E) = 


1 


The torsion scalar is also obtained to be T = 6 a' 2 /(N 2 a 2 ) = 6 P 2 , where a prime denotes 4^. 
After some algebra, action (11) reduces to 

S = M 2 J | — Na 3 (/ + T)-m 2 [/3i (3Na 2 + a 3 ) + p 2 (3iVa + 3a 2 ) + p 3 {N + 3a) ] } d A x + S m . 


(14) 

(15) 

(16) 


(17) 


For further convenience, we define some variables here, namely 


x = a~ , (18) 

2 /i {x) = 3/3ix + 3/3 2 x 2 +/3 3 x 3 , (19) 

y 2 {x) = Pi +2(/3i +j3 2 )x+ (/3 2 + /3 3 )x 2 , (20) 

2 / 3 ( 2 ) = Pi + 2(3 2 x + P 3 X 2 , (21) 

2/4(2) = -pi + {pi - 2p 2 )x + (2p 2 - p 3 )x 2 + p 3 x 3 . (22) 


Now we can get the field equations by varying action (17) with respect to N and a, respectively, 


■ 77 2 P = ^Tfr - f + T- m 2 yi{x), 


-Lp = -8 Tf TT H - T + / - 2 Tf r - 4 (1 + f T ) H + m 2 y 2 (x ), 


(23) 

(24) 


where a dot denotes fr = df/dT, frr = d 2 f/dT 2 , and p, P are the energy density and pressure of 
all perfect fluids of generic matter, respectively. Notice that we do not need to know the explicit form 
of S m in the above derivation. As is well known, the energy-momentum tensor of a perfect fluid takes a 
diagonal form in the comoving coordinates, namely 


T'4 = diag (—p, P, P, P ). 

And from the standard definition of energy-momentum tensor 

2 SS m 


rpflU _ 




(25) 


(26) 


we can easily read off p and P as 


1 SS m 
a 3 5N ’ 


P = 


1 SS m 
3Na 2 5a 


(27) 


From Eqs. (23) and (24), we can get the corresponding energy conservation equation for matter, 

p = -3H(p + P) + 3m 2 M 2 Hy 3 {x). (28) 

We can recast Eqs. (23) and (24) in the similar form as in the GR case by introducing the energy density 
and pressure of the effective DE as 


p de = M 2 [f -2 T/ t + m 2 yi{x)\ , 

P de = M 2 {(8 Pf TT + 4 f T ) H + 2 Tfr ~f~ m 2 y 2 (x )} . 


(29) 

(30) 
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Then we get the modified Friedmann equations as 

6M 2 H 2 = p total , (31) 

4 M 2 H = — ( Ptotal + Ptotal) j (32) 

in which ptotal — P + Pde, Ptotal = P + Pde- From Eqs. (29) and (30), one can obtain the corresponding 
energy conservation equation for DE as 

Pde = -3 H ( p de + Pde ) - 3 m 2 M 2 Hy 3 (x). (33) 

Combining Eqs. (28) and (33) yields 


fit ^ Ptotal') — 3fl Hptotal, (34) 

which is the usual conservation equation for total energy. We see that matter and DE are interacted 
through the graviton. We define the equation-of-state (EoS) parameter of the effective DE as 


w de 


Pde 

Pde 


4 (2T/ tt + / T ) H + m 2 y 4 (x) 
-2T f T + f + m 2 yi(x) 


(35) 


IV. OBSERVATIONAL CONSTRAINTS 

Here we are interested in the late-time behavior of f(T) massive gravity. In this section, we study two 
specific functionals of f(T) as concrete examples. These two examples will be constrained by the recent 
cosmological data. We then study the evolution of the effective EoS parameter of DE. 


A. General set-up 


Since we are interested in the late-time universe, we ignore the radiation component and only consider 
the pressureless matter, whose energy density is p m and pressure is P m = 0. Furthermore, we also 
set 02 = 03 = 0, since at the background level, the terms 3/32TO 2 a -2 and 3/3 3 m 2 a~ 3 in Eq. (29) are 
indistinguishable from the space curvature and pressureless matter terms, respectively. Having these in 
mind, and by using Eqs. (19) and (24), the effective EoS parameter of DE (35) can be further simplified to 

‘PPStt ~ /r + f/T + ^n x0 (2z/3 + 1) 

Wde (2Tf TT + /t + 1) [— 2 /t + f/T + -t^Q x q(z + 1)] ’ 


where the fractional density of the pressureless matter and the graviton are given by 


n 


m 


Pm 

6 ~MpP 


and fl x 


f3\ m 2 
2 H 2 ’ 


(37) 


respectively. Note that we always use a subscript “0” to denote the present value of corresponding 
quantity. So T 0 = T(z = 0) and D x o = flx(z = 0), with the redshift 2 defined as 2 = a -1 — 1. 
Consequently, we find the expression for the dimensionless Hubble parameter E = H/Hq, namely 


E 2 (z) = 0,0(1 + z) + n m0 (l + z) 3 + JL - 2 E 2 (z) f T . (38) 

lo 

Note that E appears in both sides of this equation, and implicitly in / and /t- 

We now briefly review the cosmological data and fitting methodology used in constraining the model 
parameters. We will perform a joint analysis of the Type la Supernovae (SNIa), the baryonic acoustic 
oscillation (BAO) and the cosmic microwave background (CMB) data to break the degeneracy between 
the model parameters. For SNIa, we use the Union2.1 dataset [77] which consists of 580 data points. 
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These data are given in terms of the distance modulus Vobs{zi)- By definition, the theoretical distance 
modulus is given by 


Vth(zi) = 5 log 10 D L {zi) + no , 


(39) 


where vo = 42.38 — 51og 10 h with h the Hubble constant Ho in units of lOOkm/s/Mpc. Here, the 
luminosity distance can be calculated as 


Dl(z) 



dz 

mpv 


(40) 


in which p denotes the model parameters. Consequently, the % 2 from 580 Union2.1 SNIa is given by 


Xsiv(p) = E 

i 


[Vobs(Zi') Vth (Zi')] 
cr 2 (Zi) 


(41) 


where a is the corresponding lcr error. Following [78], we marginalize over vo by expanding the Xsn with 
respect to vo as 

Xsiv(p) = ^ — 2 /xq.B + /ipC, (42) 


where 


Mp) = E 


[f^obs f^th 5 /^0 


o,p )] 2 


B(p) = E 


ftobsiZi) /A 0 — 0? P) 

vlobte) 


<? = E 


M) 


Eq. (42) has a minimum for vo = B/C at 


Xsat(p) = ^L(P) - ■ ( 43 ) 

Since Xsn min = Xsn min ( U P f° a constant), we can instead minimize Xsn w hich is independent of vo- 
For the data of CMB and BAO, we use the shift parameter R from CMB, and the distance parameter A 
from the measurement of the BAO peak in the distribution of SDSS luminous red galaxies, since they are 
model-independent and contain the main information of the observations of CMB and BAO, respectively 
(see e.g. [79]). The shift parameter R of CMB is defined by [80] 


R = n 


1/2 

mO 


dz 

W)' 


(44) 


where the redshift of recombination z* is determined to be 1089.90 by the Planck 2015 data [81]. On the 
other hand, the Planck 2015 data have also determined the observed value of shift parameter R 0 bs to be 
1.7382 ± 0.0088 [82], So, the x 2 for CMB is 


2 _ (R-R obs ) 2 

XCMB — ~2 

a R 


(45) 


The distance parameter A of the measurement of the BAO peak in the distribution of SDSS luminous 
red galaxies [83] is given by 


±r 

Zb Jo 


dz 


-| 2/3 


E(z) 


A = n^E(zb )- 1 / 3 


(46) 













where z b = 0.35. In [83], the value of A has been determined to be 0.469 (n s /0.98) -0 ' 35 ± 0.017. Here 
the scalar spectral index n s is taken to be 0.9741 by the Planck 2015 data [82]. And the corresponding 
y 2 for BAO is 


2 (A- A obs ) 2 

XBAO — ~2 


(47) 


So, the total y 2 is given by 

X 2 = Xsn + Xcmb + X%ao • (48) 

Then we can minimize the total y 2 to get the best-fit values of model parameters. The 68.3% and 95.4% 
confidence levels are determined by Ay 2 = y 2 — y 2 min < 3.53 and 8.02, respectively, if there are 3 free 
model parameters. 




^m() H m0 b 


FIG. 1: (Color online) The 68.3% (dashed) and 95.4% (solid) confidence level contours for the power-law f(T) 
massive gravity in the (f2 m o, b) plane (left), the (fi m o, fi^o) plane (middle) and the ( b , fi^o) plane (right). The 
best-fit parameters are indicated by the solid points. 


B. The power-law case 

In this subsection, we consider the power-law functional of f(T) first introduced by Bengochea et al. 
in [33] (see also [34]), which reads 

f(T) = aT\ (49) 

where a and b are both constants. Demanding Eq. (38) to be satisfied at redshift z = 0, we have 

« =- Y^2b - ( 6H °) ■ (50) 

Substituting Eq. (50) back into Eq. (38), the corresponding background evolution reads 

E 2 = Ha;o(l + z) + H m o(l + z) 3 + (1 — Qxo ~ Hrao) E 2b . (51) 

We note that, when Q x q = 5 = 0, this model corresponds to ACDM model in fact. There are 3 free 
parameters in this model, namely £l x o, H m o, and b. By minimizing the corresponding total y 2 in Eq. (48), 
we find the best-fit parameters U x o = 0.004, f l m0 = 0.278, and b = —0.015, while Xmin = 562.255. In 
Fig. 1, we present the corresponding 68.3% and 95.4% confidence level contours for the power-law f(T) 
massive gravity in the (D m o, b), the (D m o, H x o) and the (5, H^o) planes, respectively. We find that ACDM 
model (corresponding to = b = 0) is still consistent with the observations at the 68.3% confidence 
level. In Fig. 2, we present the evolutionary curve of Wd e with the best-fit values of model parameters. 
Apparently, the power-law case shows a phantom crossing behavior with the crossing of phantom divide 
line occurring at redshift z ~ 7.71. 
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FIG. 2: (Color online) The evolutionary curve of the effective EoS parameter of DE with the best-fit values of 
fl m o, IGo and b, for the power-law f(T) massive gravity. The left and right panels correspond to large and small 
range of redshift, respectively. The phantom divide line (Wde = —1) is also indicated by the dashed lines. 


C. The exponential case 


In this subsection, we consider the exponential functional of /(T) first introduced by Linder in [34], 
which reads 


/(T) = -aT( l-e 6 ^), 


(52) 


where a and b are both constants. Demanding Eq. (38) to be satisfied at redshift z = 0, we have 

1 D m o 

“ ~ 1 - (1 - 2 b)e b ' 

Substituting Eq. (53) back into Eq. (38), the corresponding background evolution reads 


e 2 — n x o(i + z) + n m0 (i + ~ ) 3 + 


1 9m0 D 


1 - (1 - 2 b)e b 


GO e 2 




(53) 


(54) 


We note that, when = 6 = 0, this model corresponds to ACDM model in fact. There are 3 free 
parameters in this model, namely £l x g, and b. By minimizing the corresponding total y 2 in Eq. (48), 

we find the best-fit parameters fl x o = 0.008, Q m o = 0.279, and b = 0.021, while Xmin = 562.262. In 
Fig. 3, we present the corresponding 68.3% and 95.4% confidence level contours for the exponential f(T) 
massive gravity in the (f2 m o, b), the (f2 m o, D x0 ) and the (6, fl x o) planes, respectively. We find that ACDM 
model (corresponding to = b = 0) is still consistent with the observations at the 68.3% confidence 
level. In Fig. 4, we present the evolutionary curve of Wde with the best-fit values of model parameters. 
Apparently, the exponential case shows a phantom crossing behavior with the crossing of phantom divide 
line occurring at redshift z — 1.43. 


V. CONCLUSIONS 

In this paper, we extend /(T) theory and dRGT massive gravity to a new kind of modified gravity 
model, namely f(T) non-linear massive gravity, by adding the dRGT mass term to f{T) theory. This 
mass term is formulated in the vierbein formalism to agree with the teleparallel space-time. Since the 
resulting field equations are 2nd order, and the mass term does not contain the square root of the metric 
when using vierbein formalism, this theory is easier to analyse than f(R) non-linear massive gravity. 
Besides, thanks to the rich structure of f(T) sector and massive graviton, it is natural to expect that this 
theory could also unify the early inflation and late-time acceleration in a consistent framework, and we 
leave this issue to the future works. 
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^m0 b 


FIG. 3: (Color online) The 68.3% (dashed) and 95.4% (solid) confidence level contours for the exponential f(T) 
massive gravity in the (f2 m o, b) plane (left), the (fi m o, IGo) plane (middle) and the ( b , D^o) plane (right). The 
best-fit parameters are indicated by the solid points. 


We then investigate the cosmological evolution of f(T) non-linear massive gravity. In particular, we 
study the power-law and exponential cases of f{T) massive gravity as two toy models. We then perform 
a joint constraint on the model parameters by the recent data of SNIa, CMB and BAO. We find that 
the power-law and exponential f(T) massive gravity are consistent with these cosmological observations. 
Furthermore, we explore the evolution of the effective EoS parameter of DE, and find that it can realize 
the crossing of the phantom divide line from the quintessence-like phase ( Wd e > — 1) to the phantom-like 
one (wde < —1) by using the best-fit parameters obtained from the above cosmological constraints. We 
note here that the recent data shows great possibility that the EoS parameter of DE crosses the phantom 
divide line from the quintessence-like phase to the phantom-like phase as the redshift z decreases in the 
near past [84-87]. Although there exist some complicated specific f(T) models to realize the phantom 
crossing behavior [41, 88], in general, especially in the original power-law and exponential f(T) gravity, 
phantom crossing is impossible [88-90]. So, our results are of interest. 
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